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ABSTRACT 

In the construction of multi-mass King-Michie models of globular clusters, an 
approximated central energy equipartition between stars of different masses is 
usually imposed by scaling the velocity parameter of each mass class inversely 
with the stellar mass, as if the distribution function were isothermal. In this pa- 
per, this 'isothermal approximation' (lA) has been checked and its consequences 
on the model parameters studied by a comparison with models including central 
energy equipartition correctly. It is found that, under the lA, the 'temperatures' 
of a pair of components can differ to a non-negligible amount for low concen- 
tration distributions. It is also found that, in general, this approximation leads 
to a significantly reduced mass segregation in comparison with that given un- 
der the exact energy equipartition at the centre. As a representative example, 
an isotropic 3-component model fitting a given projected surface brightness and 
line-of-sight velocity dispersion profiles is discussed. In this example, the lA 
gives a cluster envelope much more concentrated (central dimensionless poten- 
tial W = 3.3) than under the true equipartition {W = 5.9 x 10~^), as well as 
a higher logarithmic mass function slope. As a consequence, the inferred total 
mass (and then the global mass-to-light ratio) results a factor 1.4 times lower 
than the correct value and the amount of mass in heavy dark remnants is 3.3 
times smaller. Under energy equipartition, the fate of stars having a mass below 
a certain limit is to escape from the system. This limit is derived as a function 
of the mass and W of the giants and turn-off stars component. 

Key words: stellar dynamics - methods: analytical - methods: numerical - 
galaxies: kinematics and dynamics - galaxies: star clusters - globular clusters: 
general 



1 INTRODUCTION 

Since the end of 60's, models based on the Maxwellian 
phase-space distribution function (DF) have been ex- 
tensively used in the so called parametric fit of glob- 
ular clusters light profiles, in order to infer important 
properties like t he mass-to-light ratio an d the mass func- 
tion slope fsee iMevlan fc Heggid Il99^ . sect. 7.7, for a 
review). Corresponding to the isothermal distribution 
iBinnev fc Tremainelll987^ . the Maxwellian is expected 
to arise from relaxation phenomena like the interaction 
between the stars and the rapidly varying mean potential 
during the 'violent relaxation' subsequent to the clus- 
ter formation (Lvndcn-Bell 1967; Nakam ural l200(i) . or 
the close random stellar encounters (before the onset of 
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gravothermal oscillations, see iBinnev fc Tremaindll987l. 
sect. 8.4, and references therein). 

Nevertheless, the isothermal sphere has infinite mass 
and, moreover, the spatial extent of a real system cannot 
be infinite, also because of the presence of the external 
galactic field, thus some kind of cutoff for high-energy 
stars at a given tidal radius rt must be included. Among 
the vario us possibilities , the 'lowered' Maxwellian pro- 
pose d by i Michie 1* 19631) . iMichie fc BodenheimeJ lll963h 
and iKina CI 966) . on the basis of Fokker-Planck calcu- 
lations, has been the most successful one, both because 
of easy of computation and implementation and be cause 
of th e rather good agreement with observations iKind 
[l_98l'). Then, after the first results achieved with single 
component (i.e. equal star mass ) models (e.g.. IlllingworthI 
11979 : llllingworth fc Kin3ll977l) . it was clear that to re- 
produce the observed surface brightness profiles it is often 
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necessary to include stellar objects with different mass, 
i.e. multi-component DFs. 

A general form for a n-component isotropic 
DF can be given as a linear combination of sin- 
gle lowered Maxwel lians (|Da Costa fc FreemarJ llfi 
iGunn fc OriflHnllTflT^: 



/(r,«) = ^a/fe(r,t;: 
fc=i 

with 



'E/a 



1, if £ < 0, 
if £ ^ 0, 



(1) 



(2) 



being E = v'^ /2 + '^{r) the energy per unit mass and ^'(r) 
the mean total potential defined so to have *I'(rt) = 0. 
The Ck are normalization constants related with the 
global mass function, and the parameters at are con- 
nected to the velocity dispersion and depend on the stel- 
lar mass nil, 1712, ■ ■ ■ , rUn of the components (and, as we 
will see, on their concentration). 

In these models, the inclusion of the star mass as 
a further degree of freedom makes that another dy- 
namical process can be taken into account: the energy 
equipartition that is expected to occur in the denser 
regions via 2-body c ollisions between stars with differ- 
ent masses jSpitzei j Il96itl . An important consequence 
of this phenomenon, that has bee n found on various 
globu lar clust ers by means of H ST jK ing^ Sos in fc CooJ 
Il995l: ISosinl Il997j: [Ajb row. De Marchi fc Sahul 1200^ 
VLT I Andreuzzi et alTl200ll) and SDSS observations 
iKoch 6^1X1200^ is the 'mass segregation', i.e. a 
'drift' of low-mass stars towards orbits located mainly 
in the cluster outskirts while heavier stars tend to 
concentrate in the inner region. For this reason, 
these multi-mass models have been used also in A'^- 
body simulations of clusters ev olution to generate 
mass-segregated initial conditions dCornbes et alJll999l : 
ICapuzzo Dolcetta. Pi Matteo fc Miocchill2005^ . 

When accurate enough kinematical data became 
available, the so-called King-Michie (KM) models, in 
which the DF of Eq. Q includes velocity anisotropy 
through the dependence of Ck on angular moment, have 
been employed in order to reproduce the steep decrease 
in the p rojected velocity disp ersion observed in many 
clusters iGunn fc Griflinlll97gl) . Indeed, some degree of 
anisotropy in the cluster outer regions should be present 
as a reminiscence of the initial, quasi -radial, collapse 
phase of cluster formation ^SpitzeI^|l987^ . This is also an 
important issue from a stellar dynamics point of view, be- 
cause it is still unclear up to what extent this anisotropy 
can 'survive' the isotropising effect of the tidal field 
jTakahashi fc Leell200ol: iBaumgardt fc Makindl2003t) . 

The most serious disadvantage of this 'model- 
building' approach is that many, and more and more, 
degrees of freedom have to be included in order to im- 
prove t he fit with the observed profiles (see the discus- 
sion in iMerritt. Mevlan fc Mavorl Il997^ that were also 

J iroved to poo rly constraint all the various parameters 
Deionghe fc Mcrri tt 1992,'). For this reason, t he alter- 
native 'non-parametric' approach was proposed jMerritd 
Il993allbl) in which the relevant globular clusters features 



can be directly extracted by the observed profiles without 
making any assumptions on the form of the underlying 
DF^. Nevertheless, such an approach gives no direct in- 
dications on what are the actual dynamical mechanisms 
that determine the cluster evolution. For this reason, in 
our opinion, the parametric technique can still provide 
a valuable help in the comprehension of globular clus- 
ters dynamics, especially when used in conjunction with 
the non-parametric approach in order to test the validity 
of the various physical assumptions made in writing the 
form of the given DF. For instance, it still remains an 
intriguing problem that o f understanding why the fo rm 
of the giants' DF found bv lGebhardt fc Fischer! ^1995^ in 
their non-parametric study of the two non-collapsed clus- 
ters NGC 362 and NGC 3201, is clearly non-Maxwellian 
in the central region. 

By the way, it is crucial that all the dynamical mech- 
anisms taking into account in the parametric method 
are correctly and coherently 'implemented'. Neverthe- 
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just by setting the square of the velocity parameter of a 
star component as cr^ oc l/nik, which is a prescription 
that is only asymptotically correct in the limit of isother- 
mal DF. 

Various author s (e .g.. lMerritt| L98l'; 

iKondrat'ev fc Ozernovl Il982l: iMevlan fc Heggid [199711 
have already pointed out that such an 'isothermal 
approximation' does not strictly imply the central 
equipartition of kinet ic energy between stars of differ- 
ent mass. Moreover, iMerritd il98ll) drew interesting 
conclusions on the 'equipartition instability' that he 
re-considered on models with exact central equipartition. 
Nevertheless, one may still think that the isothermal 
approximation leads, after all, to a central energy 
equipartition that is, to practical p urposes, alway s 
acceptable. To this regard, following iMerritd il98H) . 
we construct multi-component models in which exact 
equipartition of kinetic energy is enforced at the centre 
(Sect. |5J and, in Sect. |3 compare their properties and 
fitting parameters with those of models constructed in 
the standard way. Discussions and conclusions can then 
be found in Sect.0] 



2 THE INCLUSION OF CENTRAL ENERGY 
EQUIPARTITION 

In Eq. 0, the non- isothermal nature of f^, due to the 
presence of the term —1, makes each component veloc- 
ity dispersion to depend on the position r, such that 
global equipartition is, in fact, impossible. By the way, 
a global equipartition would be unrealistic for globular 
clusters since in their outer regions relaxation time is sup- 
posed to be longer than their age. Nevertheless, such an 



^ Actually, some assumption must be made in order to 
infer the M/L behaviou r and the mass function form 
jGebhardt fc FischedFigOSl) . 
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equipartition is expected to occur at th e centre (r = 0) 
where the relaxation time is shortest ijSpitzer fc Hard 
ll97ll : lGunn fc GriffiJl979^ . This local equipartition can 
be formally ensured by imposing the relations: 



(3) 



where the first mass class is chosen as the 'reference' com- 
ponent, without the loss of generality, and denotes the 
modulus square of the velocity of the stars in the A;-th 
component at r = 0, i.e. 



{■"k) = / fk(0,v)v'^dv 



fk{0,v)v dv 



(4) 



'0 /JO 

with w = a/2|*o| the escape velocity, being 'I'o = ^'(0) 
the central potential. 

In the following, Wk = I'Pol/crfc denotes the dimen- 
sionless central potential and, moreover, when referring 
to a quantity of the generic k-th component, for simplic- 
ity the index k is omitted (if unambiguous). 
It can be shown that: 



l)v'^dv 







2W 



(e e ' 



l)v'^dv 



^e^Vy-^/^erf(VW/) ■ 



4W2 



(5) 



with erf (x) = (2/0F) T e"* dt the error function. Using 
the relations (|KJ and Q, Eq.Q can be rewritten, after 
some terms rearrangement, as 



(v^) = (j'^Wk[W), 
where one defines 
k{W) 



(7) 



— — e 



erf(^l^) - ^W^^^ 



W 



3/2 



5 



3 



(8) 



Hence the constraint imposes that ma^Wn{W) — 
mialWiK,{Wi). Therefore, being alWi = a'^W = l^o], 
the central equipartition leads to the following relation 
for W: 



mK{W) = miK(Wi). 



(9) 



Given Wi and m/mi, the quantity W can be found for 
any component provided the knowledge of the inverse 
function of k, which exists for W > and < k < 
6/7 - being k{W) a monotonically decreasing continuous 
function (Fig. - and which can be easily computed 
by numerical means (e.g. using the bisection algorithm). 
Thus, 

W = ( — k{Wi)] . (10) 
Km J 

Then, once given ai, the parameter a for any other com- 
ponent is determined by the relation: 




Figure 1. Plot of the function k{W') defined in Eq.JSJ. Note 
tliat k(W) ^ 6/7 for W ^ 0+ . 



2W1 



(11) 



Note that, since k < 6/7, Eq. @ implies the existence of 
a lower limit for the stellar mass in the components, in 
fact one has that 

m > ZK(W^i)mi. (12) 

This is just a direct consequence of the obvious condition 
{v^) < and it means that, if exact equipartition is 
achieved in the central region of a real cluster, objects 
with a mass lower than the limit imposed by condition 
acquire enough energy to escape from the system. 
Hence, they can be present there only if they are bound in 
multiple systems formed before 2-body relaxation takes 
place (eg., light remnants in primordial binaries). In this 
respect, note that for a component with such light ob- 
jects, it might be appropriate to multiply the fk of Eq.Q 
by a factor that depopulates the phase-space region with 
low angular momentum L, so as to make its inclusion in 
the model consistent with the equipartition. A suitable 
factor could be, e.g.. 



exp (-LV2rLai) 



(13) 



being the 'avoidance radius', rav, a further free param- 
eter. Notice, finally, that with the typical values mi = 
0.7 Mq and Wi = 8, Eq. lO yields m > 0.3 M©. 



2.1 The isothermal approximation 

As known, the central dimensionless potential W de- 
termines the 'form', the shape, of the density distribu- 
tion given by the fk of Eq. being the other param- 
eters (e.g. a and p(0)) just scaling factors. In particu- 
lar, W is directly related to the King concentration co- 
efficient c — log(rt/rK), with tk = 3a{4nGpo)~^^^ the 
King radius; they satisfy the approximate linear relation 
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c ~ 0.06 + 0.2W jKindll966l : lBinnev fc Tremainj|l98'i1) . 

Then, when ^ 1 it is clear that the presence of the 
tidal cut will have a negligible role, in fact one can say 
that fk — exp{—E/al) at the centre and for small veloci- 
ties {v^ <C Wa^), hence the approaches an isothermal 
DF whose velocity dispersion is {v'^) — 3a^. Thus, the 
'isothermal approximation' (hereafter lA) yields, from 
Eq. (|3J, the commonly adopted equipartition condition 



0-1 • 

m 



(14) 



which allows to determine cr for any component and, in 
turn, to evaluate its central dimensionless potential as 
WiA = Wiaf/a'^ = mWi/mi. 

Note that under the lA, Eq.JSJ correctly leads to the 
same condition since the term dominates for ^ 1, 
giving k{W) ~ 3/W, hence Eq. @, for W > 1 and 
Wi S> 1, gives just W ~ mWi/mi. 

It is now worth checking and quantifying more pre- 
cisely the goodness of the I A: let us evaluate the actual 
'temperature' ratio T/Ti between a generic component 
and the reference one at the centre of the cluster, and 
compare it with the unitary value actually demanded by 
the equipartition. From Eq. Q we have that 



mK,{Wij 



mK{rnW\/m\) 



(15) 



(16) 



mi{vf) mi«:(Wi) miK,{Wi) 

In Fig|5| this ratio is plotted for various values of Wi 
as a function of the mass ratio m/mi. As expected, the 
isothermal approximation works better, i.e. T/Ti is closer 
to the unit, for higher Wi and for m > mi, because in 
these regimes the DF of both components are closer to 
the isothermal distribution. However, although one can 
affirm that the dynamical condition of central 'thermal 
equilibrium' is approximately valid within a limited range 
of masses and for sufficiently concentrated models, nev- 
ertheless the model structure is rather sensitive to the 
temperature ratio. Let us consider the quantity 

WiA _ mWi 
'W ^ miW 

where W is given by Eq. IIOII . This ratio quantifies the 
change affecting the fundamental structural parameter of 
the clusters profile, when the lA is adopted instead of the 
exact energy equipartition. 

As shown by Fig|3 one can see that Wia can be very 
different from W , especially for low concentration models, 
e.g. for non-collapsed globular clusters {Wi < 10). In 
particular, given a pair of components, the use of the lA 
implies a higher concentration for the component with 
lighter stars, and a lower concentration for the heavier 
stars. In other words, with respect to the correct energy 
equipartition assumption, the lA leads to a reduced mass 
segregation. 

Then, as a consequence on the parametric fitting of 
globular clusters profiles, one expects that, once adjusted 
the distribution of the brightest stars component (giants 
and turn-off stars) to fit the observed luminosity profile, 
the lA gives: (i) a steeper and too much concentrated 
mass density profile in the outer regions, being this deter- 
mined by the lighter (dwarf main sequence) stars mostly; 
(ii) a flatter density profile for the massive remnants in 
the core region; (iii) a steeper radial velocity dispersion 




I I I I I I I I I I 

-0.4 -0.2 0.2 

log(m/mJ 



0.4 



Figure 2. Temperature ratio at the cluster centre between a 
generic component with mass m and a 'reference' one having 
mass mi and Wi as labelled, when the energy equipartition is 
imposed using the 'isothermal approximation'. The thick Hne 
marks the low-mass cutoff given by 1121 . 



profile for the giants, as a consequence of the steeper de- 
creasing of the potential. Point (i) implies then a lower 
total cluster mass, and so a lower mass-to-luminosity ra- 
tio, in respect to the exact energy equipartition. 

Notice that all this discussion remains valid also in 
the case of KM models where the presence of velocity 
anisotropy is considered. Indeed, in these models, the co- 
efficients Cfe in Eq.(0 are equal to at exp{—r^ v'j^ /2r^al) , 
where Ta is the length-scale above which anisotropy is im- 
portant and v± is the further state variable corresponding 
to the modulus of the tangential velocity (see Appendix 

Thus, for r = there is no anisotropy and C'k = O-k 
is independent of u^. 



3 CONSEQUENCES ON THE OBSERVABLE 
PROFILES AND RELEVANT 
PARAMETERS 

To understand better the effects of the lA on the mod- 
els morphogical and kinematical features, a simple case 
of an isotropic (ra = oo) 3-component model is consid- 
ered. Though representing a crude simplification in re- 
spect to real clusters with a continuous mass spectrum, 
it is sufficient to provide significant indications on the 
consequences of the approximation we are discussing. 
The model is made up of: (i) giant and turn-off stars 
(mi = 0.75 M0, total mass Mi), determining the light 
profile; (ii) main sequence dwarf stars and light remnants 
(m2 = 0.26 Mq, M2/M1 = 22), considered to be dark, 
but dominating the gravitational potential; (iii) heavy 
dark remnants (ma = 1.2 M©, M^/M-i — 0.28, remnants 
of objects with mass in the range 4-8 M0 ). All is consis- 
tent with a mass function dN oc m~'^d log m with x = 1.7, 
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Figure 3. Ratio between the dimensionless central potential 
evaluated under the lA of a component with mass m, and 
that turned out by the correct energy equipartition with a 
'reference' component having mass mi and various Wi. The 
vertical asymptotes correspond to the m lower limit 1121 . 



and with condition I12I I. The surface brightness profiles 
correspond to a concentration parameter c ~ 1. In Fig.|l] 
the surface brightness and the projected line-of-sight ve- 
locity dispersion profiles given adopting the lA and the 
exact equipartition, are shown for a comparison. The ra- 
dius is given in unit of the core radius, Vc, at which the 
luminosity drops to 1/2 the central value. The parameter 
Wi has been adjusted in such a way to give similar light 
profiles in both cases. 

Thus, confirming the points discussed at the end of 
Sect. 12.11 the lA gives a cTp profile appreciably steeper 
than in the correct case. The cluster envelope is also 
too concentrated {W2 = 2.3 instead of the correct W2 = 
8.7 X 10~^) and the massive remnants have a much flat- 
ter density profile in the core region {W3 — 11 instead 
of W3 = 16). Though to a less extent, this occurs also 
in the presence of velocity anisotropy, as shown in Fig|^ 
Moreover, the difference is lower also when a less peaked 
mass function is adopted, as indicated in Fig. |S] that 
refers to the same isotropic model above-described, but 
with X = 1.2, i.e. with: mi = 0.75 M©; m2 = 0.28 Mq, 
M2/M1 = 14; ms = 1.2 Mo, M3/M1 = 0.76. Of course, 
in this case the reduced discrepancy is due to the less 
importance of the lightest component in determining the 
overall potential. 

Now, to quantify the change affecting all the relevant 
parameters because of the use of the lA, these have to be 
compared with the fitting parameters that, under the cor- 
rect energy equipartition, give the same projected surface 
brightness and velocity dispersion profiles. In Fig. |7|the 
profiles corresponding to a typical non-collapsed cluster 
(c = 1.2 in this case), are plotted. The resulting param- 
eters are written in Table Besides the difference in 




Figure 4. Upper panel: projected velocity profile for the lumi- 
nous component of the isotropic 3-component model described 
in the text. Lower panel: projected surface brightness profile. 
Solid line: profiles obtained in the case of exact energy equipar- 
tition among the components; dotted line: profiles obtained 
under the I A. 

the lightest component concentration, one can immedi- 
ately see the discrepancy in the total mass (and so in 
the mass-to-light ratio) and in the mass function index. 
In this example, the lA leads to a 1.4 times lower total 
mass and to a steeper mass function, x = 1.5 instead of 
1.2, giving an amount of mass in heavy remnants which is 
3.3 times lower. See AppendixIXlfor a detailed description 
of the procedure followed to construct the KM models. 



4 CONCLUSION AND DISCUSSION 

In the multi-mass King-Michie models usually employed 
for the study of morphology and kinematics of globular 
clusters, the energy equipartition at the cluster centre, 
between stars of different masses, is imposed by setting 
the square of the velocity parameter of the fe-th mass 
class, o"! in the DF of Eq. Q, inver sely proportion al to 
the mass, ruk, of that compon ent. As lMerrifrg lll98ll) and 
iKondrat'ev fc Ozernovl il982f) emphasized, this is correct 
only in the limit of isothermal distributions (i.e. with 
central dimensionless potential W — > 00). This 'isother- 
mal approximation' (lA) has been checked and the conse- 
quences studied in this paper. To this purpose, a general 
procedure to construct n-component models with exact 
energy equipartition at the centre is used and described 
in detail. 

From the point of view of the 'thermal equilibrium'. 
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Table 1. Fitting and derived parameters for the 3-component isotropic model whose profiles are plotted 
in Fig. |7| The reported parameters are: the mass function logarithmic slope (x) and the system total mass 
(M, in arbitrary units); then, for the components: the total mass (M^, in arbitrary units), the central 
dimensionless potential (Wk), the velocity parameter {cr^., with ai = 1) and the temperature ratio. In both 
cases we have mi = 0.75, m2 = 0.28 and ms = 1.2 Mq. 
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Figure 5. The same as in Fig.l^but with an anisotropy radius 
fa = 4rc 

it is found that the 'temperatures' of a pair of compo- 
nents can differ to a non-negligible amount for low con- 
centration distributions. For instance, taking a giant and 
turn-off stars component with mass mi ~ 0.8 M© having 
a concentration c ~ 1, along with a typical dwarf main- 
sequence star component with ma ~ 0.2, one has that, 
under the lA, this latter component has a temperature 
(m{ii^) /2) which is actually a factor 0.75 lower than that 
of the giants. 

To investigate the consequences on the observable 
profiles, a simple 3-component isotropic model was gen- 
erated resembling roughly that of a typical non-collapsed 
cluster. Aside from the two above-mentioned compo- 
nents, heavy remnants with ma = 1.2 M© were also in- 
cluded, while the components total mass were varied ac- 
cording to different mass function logarithmic index x. 
Once adjusted the giant profile to match the same hght 
distribution, it is found that (with x — 1.7) the cluster 
envelope is much less concentrated in the exact equiparti- 
tion case {W = 8.7 x 10"^) than under the lA {W = 2.3). 




Figure 6. The same as in Fig. |l]but with a flatter mass func- 
tion (x = 1.2). 

Moreover, the lA gives a flatter density profile for the 
massive remnants in the core region {W = 11 instead of 
W = 16). As regards the velocity dispersion profile, the 
approximation leads to a steeper profile for the giants, as 
a consequence of the steeper decreasing of the potential. 
However, it is found that both the presence of anisotropy 
and lower mass function slopes reduce the variations in 
respect with the correct case. 

Then, in order to evaluate the change in the most rel- 
evant parameters, these were fixed so as to give the same 
projected surface brightness profile and the same line-of- 
sight velocity dispersion curve, with models constructed 
both with the correct energy equipartition and employing 
the lA. It is found that the inferred total mass (and, cor- 
respondingly, the global mass-to-light ratio) results 1.4 
times lower than the correct value, with a steeper mass 
function that leads to a 3.3 times smaller amount of mass 
in heavy dark remnants. 

By the results of this work it can be stated that the 
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Figure 7. Projected velocity dispersion (upper panel) and 
surface brightness (lower panel) profiles for the luminous com- 
ponent of the isotropic 3-component model made under the lA 
(dotted line) and with the correct energy equipartition (solid 
line) . They give similar curves with different fitting parameters 
(see Table 



lA leads to multi-mass parametric fit of non-collapsed 
globular clusters profiles, corresponding to models not 
consistent with tiie assumption of energy equipartition 
at the centre. Consequently, the estimate of many impor- 
tant parameters (such as: total mass, mass-to-light ratio, 
mass function slope, heavy remnant abundances, etc.) is 
appreciably affected by the non-unitary ratio among cen- 
tral components temperatures. Moreover, in the multi- 
mass models used so far, condition (I12II is usually ig- 
nored. Indeed, it is assumed the presence of centrally 
concentrated light objects (hght remnants or dwarf main 
sequence stars, with mass as low as 0.1 M©) that, on the 
contrary, should not be present in the central region, since 
they acquire there sufflcient kinetic energy to leave the 
system. 

Thus, though the models with lA can reproduce 
rather well the observable profiles, they are inconsistent 
with one of the dynamical processes (energy equiparti- 
tion) they rely on. This is incoherent with the underly- 
ing 'philosophy' of the parametric approach. By the way, 
if some dynamical constraints were relaxed without any 
theoretical justification, then non-parametric techniques 
would definitively be more appropriate. 

Whether or not the energy equipartition takes ac- 
tually place at the cluster centre is an entirely differ- 
ent problem (to this regard, our opinion con trasts to 
what claimed in iKondrat'ev fc Qzernovlll982D . framed 
into a still controversial subject. In this respect, it is 



worth noting that IPrvor et alj il99ll) deduced, thanks 
to a parametric study of low-concentration clusters, 
a relatively high and sharp cutoff in the mass func- 
tion, at 0.3-0.4 M0, for NGC 5466 and M55, which 
is in goo d agreement with condition 11211. Fur ther- 
more, th e [P^^sce^De Marchi fc Romaniellol lHoOsI) and 
iDe Marchi et al.1 (Il995l) HST deep photometric studies, 
concerning with the region around the half-mass radius 
of various clusters, indicates a general remarkable flat- 
tening of the mass function below 0.3 Mq. This was 
also observed, for higher masses, in the core of 47 Tuc 
jParesce. De Marchi fc Jedrzeiewskill995l) . Then, energy 
equipartition is, at least, not contradicted by these obser- 
vational evidences about the lacking of low-mass stars in 
the central regions of clusters. 

From a theoretical point of view, besides the effi- 
ciency of 2-body collisions, also the relative importance 
that the different components have in determining the 
total gravitational potential influences much the equal 
distribution of kinetic energ y among them and, as flrstly 
envisaged bv ISpitzeil 1^6^, it could be impossible, un- 
der certain conditions, to be attained. In a few words, 
this happens when a much more massive stellar species 
(with mass 7712) has a negligible influence on the overall 
potential determined by a much lighter one (mi ^ 7712, 
but Ml ^ M2), so that the heavier stars get more and 
more concentrated (while the others expand), becoming 
a self-gravitating system with a negative heat capacity. 
Thus, the more kinetic energy it loses, the higher will be 
its temperature departing faster and faster and more and 
more away from the energy equipartition with the lighter 
component, leadi ng to the so-calle d 'equ ipartition in- 
stability'. Though llnaeaki fc Saslawl 1^83) showed that 
equipartition never oc curs in their Fokker-Plank models, 
ISpitzer fc Hart ||l97l|) and the more recent M onte Carlo 
simulations by Walters. Joshi fc Rasid ll200Cl) have basi- 
cally conflrmed the original Spitzer's criterion for a 2- 
component model. In this respect, one must note that 
neither the relation 



/MaX /ma 
\mJ Vmi 



s; 0.32 



(17) 



suggested by IWatters et al.l i2000l) for equipartition to 
occur, nor the Spitzer's one, are satisfied for the 2nd and 
the 3rd components of our e x ample models. We must 
remind, however, the iMerritti il98 J) conclusion, drawn 
by re-considering the Spitzer's argument on models with 
exact central equipartition, affirming that this latter can 
be reached even if Spitzer's condition is violated. Perhaps, 
more convincing indications could be provided, in the 
near future, by direct and collisional A''-body simulations 
with a sufficiently large A'^ and, more important, with a 
realistic mass function ( see the introductory discussion in 
llnaeaki fc Saslaw^ll985^ . 

To conclude, it can be stated that in the parametric 
study of low-concentration clusters, it appears appropri- 
ate to take particular care to the inclusion of the energy 
equipartition process in KM multi-mass models, since the 
structure of each star component and, consequently, all 
the relevant parameters, are significantly influenced by 
the presence of thermal equilibrium at the centre. Fur- 
thermore, the coherence of the parametric approach must 
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be preserved to make significant comparison with the re- 
suhs of non-parametric (or other) techniques. 
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APPENDIX A: THE CONSTRUCTION OF THE KING-MICHIE MODEL 

Let us consider the anisotropic King-Michie DF 

n n 

f{r,v,rvj_) ^^fk{r,v,rv^) = ^ Qfe exp(-r^t)i/2r^(jfe)/fe(r, -y), (Al) 
fc=i fc=i 

where fk is that of Eq.Q and v± = usinry, with rj the angle between the velocity and the radial direction. To construct 
a self-consistent model it is necessary that the potential \l/(r-) in the DF is that generated, according to the Poisson's 
equation 

^ + -r^-^-Gpir), (A2) 
by the mass density 



pir) = ^ pk = 2tt ^ mfe 

k=l fc=l 



fkiv /2 + 'i{r),rv sin rj) sin r]dT] 



v^dv (A3) 



determinated, in turn, by the DF itself. As done in lOunn fc GriffinI l)l97!j) . the problem is solved through the numerical 
integration of Eq. IIA2II . with the initial conditions ^'(O) = —Wiai, *'(0) = 0, employing a 4-th order Runge-Kutta 
method. However, a completely analytical expression was used for the density, without the needing of numerical quadra- 
tures. Indeed, it can be shown that, for the fc-th component, the integral in Eq. I|A3|I gives: 



Pfe(3',r) = ak : 



.(A4) 



where Wfc = — 'I'/cr^, f = r/ra and erfi(a;) = (2/^7r) Jg" e* dt is the imaginary error function that can be computed 
through truncated power series expansion, so as to have: 

r ^1 ^ / 2e--'S^_Jj!(2j + l)]-ia;2^-+i, if a; ^ 2, 

Vvre ernia;) ~ < ^j-u'- (A5) 

I ^;^„(2j-l)!!2-^:c-(2^+i), ifa;>2 ^ ^ 

that, with p = 4, gives a truncation error < 1 per cent for x ^ Q. On the other hand, being erfi(a;) ~ 
{2/y/^){x - 2x^3) for a; ^ 0, it is easy to verify that, for f <^ 1, Eq. I|A4|I tends to the correct expression 
for the density in the isotropic case ijBinnev fc Tremainel Il98^ . eq. 4-131). Moreover, the asymptotic behaviour 
Pfc('I',r) ~ 87rafemfc(%/2crfc)''W^^''^/15, for <^ 1, can be used to avoid round-off error in numerical evaluations. 
The input needed to construct a complete self-consistent multi-mass model with anisotropic velocity is, then: 

the components star mass {mi, m2, . . . , rrin}; 

the components total mass {Mi, M2, . . . , M„} and the mass-to-light ratio {Ti, T2, . . . , T„}; 
the velocity parameter of the 'reference' component (ti; 
the central density of the ref. component po.ii 

the dimensionless central potential of the ref. component Wi (which gives 'I'o = — Wicri); 
the anisotropy radius ra (assumed to be the same for all components). 

Then, in order to impose energy equipartition at the centre, (Jk for the other components has to be fixed by the 
relation ak = aiWi/W, being W given by Ea. llUl with m — m^. On the contrary, if one wants to test the isothermal 
approximation then ak = miai/mk- 

The first normalization constant ai is readily fixed by the requirement pi(*o,0) = po,i, while the others are 
constraint by the given total component masses. Since Mk depends, through pk, on the overall potential which, in turn, 
depends on all the normalization constants, these contraints correspond to a system of n — 1 integral equations for the 
set of akyi- As usual, this system is solved iteratively, starting from an initial guessed set of ak, e.g. ak = aiMk/Mi, 
halting the process when the difference between two subsequent sets of ak is within a given tolerance. 

At the end, given p{r) and the potential ^'(r), all the relevant quantities can be evaluated using also Eq. Ij Alll. 
Thus, the projected line-of-sight velocity dispersion of the fc-th component is given by (iRichstone fc TVemainelll984l) : 



— +7^ 



mkfk{r,v,rv±)v'^dv (A6) 



where x = cos 77, v± = v{l — x^)^''^ and Sk{R), the projected surface density, is 



SdR)-2l p^^^dr. (A7) 



Of course, the projected surface brightness will be given by I{R) = X/fe-i '^k ^Sk{R)- 
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